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Can AGN feedback-driven star formation explain the size 
evolution of massive galaxies? 
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ABSTRACT 

Observations indicate that massive galaxies at z ~ 2 are more compact than galaxies 
of comparable mass at z ~ 0, with effective radii evolving by a factor of ~ 3 — 5. This 
implies that galaxies grow significantly in size but relatively little in mass over the 
past ~10 Gyr. Two main physical models have been proposed in order to explain the 
observed evolution of massive galaxies: 'mergers' and 'puffing-up' scenarios. Here we 
introduce another possibility, and discuss the potential role of the central active galac- 
tic nucleus (AGN) feedback on the evolution of its host galaxy. We consider triggering 
of star formation, due to AGN feedback, with radiation pressure on dusty gas as the 
driving feedback mechanism. In this picture, stars are formed in the feedback-driven 
outflow at increasingly larger radii and build up the outer regions of the host galaxy. 
The resulting increase in size and stellar mass can be compared with the observed 
growth of massive galaxies. Star formation in the host galaxy is likely obscured due to 
dust extinction and reddening. We suggest a number of observational predictions of 
our model, and discuss possible implications for AGN feedback-driven star formation. 
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1 INTRODUCTION 

Recent observations confirm that massive (M* ~ 10^^ A/©) 
quiescent galaxies at high redshift {z ~ 2) are more 
compact than galaxies of comparabl e mass a.t z ~ 
JPaddi et al]|2005l: iTruiillo et al ll2007l: Ivan Dokkum et all 



20081: iBezanson et al.l l2009l : Ivan Dokkum et al.l \20m . 



Patel et al 



20121 ). At a given stellar mass, the efi^ective radii 
of compact, early- type galaxies at 2: ~ 2 can be a factor of 
~ 5 smaller than their c ounterparts in the local universe 
(|van Dokkum et al.ll200^ ). While massive compact galaxies 
seem to have been common in th e past, such dense o bjects 
are found to be very rare today (|Truiillo et al.l [20091 ). This 
implies that massive galaxies have undergone a substantial 
evolution, with a strong increase in size, over the last ~10 
Gyr. Although a certain spread in the sizes of compac t 
galaxies is observed at high redshifts jSzomoru et al.ll2012l ). 
and slight biases against massive a nd compact object s at 
low redshifts cannot be excluded (|Tavlor et all I2OI0I ). a 
significant evolution in the galaxy properties is still required 
between « ~ 2 and the present. 

In particular, the central densities of the high redshift 
compact objects are found to be comparable to those of 
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present-day massive ellipticals, sug gesting that the main 
growth takes place at outer radii ([Bezanson et al. I I2OO9I : 
Ivan Dokkum et al.ll201ol : |Patel et al.ll2012l ). According to the 
'inside-out' growth scenario, the central core of a galaxy is 
formed early on, while the subsequent growth occurs in the 
outer regions, with a gradual build-up of an extended en- 
velope. A remarkable characteristic of the evolution of mas- 
sive galaxies is that galaxies grow significantly in size but 
much less in mass. Indeed the evolution in effective radius 
is observed to be much stronger than the evolution in stel- 
lar mass. The combination of strong increase in size and 
little increase in mass provide important constraints on the 
possible physical mechanisms at the origin of the observed 
evolution. 

A number of physical models have been proposed in 
the literature in order to account for the evolution of massive 
galaxies from z ~ 2 to the present. The two main models dis- 
cussed so far are the 'merger' and the 'puffing-up' scenarios. 
Major mergers involve mergers of galaxies with comparable 
mass, which lead to a linear increase in both size and mass. 
In contrast, in minor mergers the primary galaxy accretes 
lower mass secondary companions, and the increase in size 
results to b e more efficient co mpared with the case of ma- 
jor mergers l|Naab et al.ll2009l ). In the puffing-up model, the 
size growth is obtained via a rapid mass loss from the cen- 
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tral r egions, attributed to quasar feedback (|Fan et al.ll2008l . 
|2010| ). 

It is now widely agreed that supermassive black 
holes reside at the centre of all massive galaxies. 
Direct relationships between the characteristic proper- 
ties of the central black hole and the host galaxy 
have been observed, including the correlation between 
black hole mass and stellar velocity dispersion and/or 
bulge mass fKor mendv fc Richsto ne 1995, 'M agorrian et al] 
19981: ICebhardt et al I l2000l: iFerrarese fc MerrittI bOOoF 



Marconi fc Huntll2003l : iHaring fc Rixll2004 ). These empir 



ical correlations are generally interprete d in t erms of the 
active galact i c nucleus (AGN) feedback ijSilk fc Rccs 19_9^ ; 
lFabianlll999l : lKin!dl2003l : iMurrav et aklbOOSl : IFabianl 12013 . 
for a review), which is expected to play a major role in the 
evolution of the host galaxy. Here we discuss the possibility 
of the AGN feedback triggering star formation in the host 
galaxy, a s cenario which we hav e briefly explored in a previ- 
ous paper l|lshibashi fc Fabianll2012l . hereafter Paper I). As 
all massive galaxies harbour massive black holes at the cen- 
tre, they should have experienced some active phases during 
their lifetime, due to accretion onto the central black hole 
and associated feedback. However, given the actual AGN 
duty cycles, most galaxies in the observational samples will 
be detected in quiescent phases. In the following, we consider 
the active state, and try to interpret the observed evolution 
of massive galaxies within the feedback-triggered star forma- 
tion scenario. In this paper, we explore the main concepts 
of the AGN feedback-driven star formation model, while the 
details of the model will be developed later and presented 
in a future work. 



2 FEEDBACK-TRIGGERED STAR 
FORMATION 

We briefly summarise here the main features of our model, 
which is discussed in more detail in Paper I. We consider 
the effects of AGN feedback on the surrounding interstel- 
lar medium of the host galaxy, and in particular the direct 
effects of radiation pressure on dusty gas. Radiation from 
the central black hole is absorbed by dust grains embedded 
in the gas, and the radiative momentum is transferred to 
the gas through the efficient coupling between gas and dust 
components. We assume that radiation pressure on dusty 
gas sweeps up ambient material into a shell, and we follow 
the evolution of the outflowing shell in the gravitational po- 
tential of the galaxy. The equation of motion of the shell is 
given by 



(1) 



where L is the luminosity of the central source, Mg{r) 
the enclosed gas mass, and MoKiir) the dark matter 
mass. In Paper I, we analysed the exact conditions for 
the escape or trapping of the shell in the galactic halo 
for different underlying dark m atter potentials (see also 
iMcQuillin .^McLaughlinl (|2012l )). In general, once the cen- 
tral black hole reaches a critical mass, or correspondingly 
a critical luminosity, the shell is able to expand outwards 
indefinitely. 



In the case of an isothermal potential, the expanding 
shell reaches a constant asymptotic velocity at large radii, 
which is independent of the initial conditions: 



GL 

2/,ca2 



2a' 



(2) 



where a is the velocity dispersion and /g the gas mass frac- 
tion where the gas mass is assumed to be a fraction /g of the 
dark matter mass. We may take as fiducial numerical values 
a typical velocity dispersion of a ^220 km s"'^, a plausible 
gas mass fraction of /g ~0.1, and central luminosities of a 
few L ~ 10"^ erg s'^ 

The passage of the outflowing shell induces squeezing 
of the surrounding inhomogeneous interstellar medium. The 
resulting compression of the cold gas may drive gas clouds 
beyond their Jeans limit, inducing cloud collapse. Thus local 
density enhancements, caused by the passage of the outflow, 
may in turn lead to a triggering of star formation. We note 
that star formation may only occur in regions of enhanced 
density (dense clumps) within the outflow, which may have 
lower speeds compared to the mean shell velocity. The star 
formation rate, i.e. the global rate at which gas is turned 
into stars in the expanding shell, is parametrized as 



tfiow(r) 



2/ga^ 



v(r) . 



(3) 



where is the star formation efficiency, and tfiowij') = 
is the local flow time. The observed star formation efficiency 
is of the order of a few percent, typically in the range ^ 
0.01 — 0.1. The second equality is obtained in the case of 
the isothermal potential in which the gas mass is directly 
proportional to the radius. We see that the star formation 
rate scales with the shell velocity, and is constant at large 

radn: M*,oo = M*(r ^ oo) ~ ^^v^. 

As the shell expands outwards, new stars are formed 
at increasingly larger radii. As long as their initial velocity 
does not exceed the local escape velocity, newly formed stars 
will remain bound to the galaxy, and follow closed orbits in 
the galactic potential. We see that feedback-driven star for- 
mation mainly adds stars in the outer regions of the galaxy, 
thus contributing to the outer growth of the host. 



3 SIZE AND MASS GROWTH 

New stars are successively formed in the outflowing shell 
according to the star formation rate given in Eq. The 
resulting stellar density distribution is approximately given 
by Mi,{r) ~ J '^ai^) jj^ simple approach, assuming 
that the gas density distribution can be parametrised by a 
power law, pg{r) = Ar~' , the gas mass distribution follows 
^ai''') ~ 47r J pg{r)r^dr oc r'^'^^. The corresponding stellar 



mass profile is given by M^(r) oc r''^ , i.e. the stellar dis- 
tribution follows the underlying gas distribution with the 
same power law slope. In particular, if the gas follows an 
isothermal distribution, Pg{r) cx r~^, the resulting stellar 
distribution will also hold the characteristic isothermal pro- 
file, p*(r) cx r^^ . A recent study suggests that the total den- 
sity profiles of early-type galaxies follows a nearly isothermal 
profile, p(r) oc (with 7 ~ —2), and that t his particular 
slope can be regarded as a natural attractor (|Remus et al.l 
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I2012I ). Detailed future work is required to see if this applies 
to the model presented here. 

Once formed, stars decouple and eventually drop out 
from the outflowing shell. The driving term due to radiation 
pressure vanishes in the equation of motion (Eq.[T]), and we 
are left with only the gravitational term. Thus the newly 
formed stars start to follow their own ballistic orbits in the 
gravitational potential of the galaxy. For simplicity, we may 
assume that stars initially share the properties of the ex- 
panding shell in which they are formed, with the initial ve- 
locity corresponding to the local shell velocity at the forma- 
tion radius. The initial conditions are given by: r*,o = '"shoii 
and u*,o = «shcii(''*,o)- In a first stage, stars continue to 
move outwards following their initial radial motion. But un- 
der the effects of gravity, they will soon be slowed down, and 
later fall back, accelerating towards the centre. As a result, 
stars orbit back and forth in the gravitational potential of 
the galaxy, and on average spend most of their time where 
their velocity is lower. For simplicity, we may also assume 
that the star's trajectory is symmetric with respect to the 
centre. We can then define the mean or average radius of 
the orbiting star as 



Ir{t)dt 
fdt 



(4) 



In the isothermal potential, the resulting average radius can 
be written as a multiple of the initial radius: f = K x r*,o, 
where the normalisation constant K depends on the initial 
velocity. For initial velocities in the range a < Voo < '2o-, we 
obtain that the corresponding mean radius lies in the range 
r*,o < r < 2r*,o- Evidently, a high initial velocity leads to a 
larger value of K, and hence a larger mean radius, up to a 
doubling of the initial radius. Further, a high initial velocity 
also implies that the expanding shell had already travelled 
to greater distances. In fact the maximal radius reached by 
the outflowing shell within a given AGN activity timescale 
is roughly given by: rmax =ro + v^- /SIagn ~ ^"00 • AtAON, 
where ro <^ Vao ■ AtAON- Thus a high shell velocity implies 
both a large initial radius and large initial velocity, which 
couple to give a large mean radius for the orbit of the star. 
The formation of new stars in the outer regions can therefore 
lead to a substantial increase in the galaxy size. 

In physical terms, the increase in shell velocity, which 
ultimately leads to the increase in galaxy size, can be 
linked to the decrease of the gas mass fraction over cos- 
mic time. Indeed high-redshift {z ~ 2) galaxies are observed 
to be much more gas-rich compared to present-day galaxies. 
Measurements of the baryonic gas fraction {fg, defined as 
the ratio of gas mass to the sum of gas and stellar masses), 
yield values of the order of ~ 0.5 in normal star-forming 
gal axies, close to the peak epoch of galaxy asse mbly activ- 
ity ijPaddi et al.ll2010l : iTacconi et al.ll2oTol . l2012l ). These ob- 
served gas fractions are at least a few times higher than 
the typical values measured in local counterparts, suggest- 
ing a significant evolution in the gas mass fraction. The 
gas phase thus accounts for a major fraction of the bary- 
onic mass in high-redshift galaxies, with gas masses being 
comparable to stellar masses. The large molecular gas con- 
tent of these gas-dominated galaxies has been suggested as 
a possible explanation for the higher specific star forma - 
tion rates observed at high redshifts (|Tacconi et al.l[201Cf ). 
Interestingly, new Hubble Space Telescope (HST) observa- 



tions reveal the presence of equally significant amounts of 
cool gas in the haloes of early- type galaxies (| Thorn et al.l 
l2012h . The detection of cold, bound gas in early-type galax- 
ies seems to be comparable to that of star-forming galaxies. 
This implies that passive galaxies also contain important 
reservoirs of cold circumgalactic gas, which can potentially 
be re-accreted and eventually made available for star for- 
mation. Therefore the quiescence of 2: ~ 2 galaxies, and 
associated quenching of star formation, may not be entirely 
attributed to the lack of cold gas reservoir. 

A rough order of magnitude estimate of the total stellar 
mass added by the feedback-driven star formation process 
can be obtained by integrating the star formation rate (Eq. 
[S]) over the AGN feedback timescale. Assuming radiatively 
efficient accretion, we may take the Salpeter time as the 
characteristic feedback timescale. We obtain typical star for- 
mation rates of >100 Mq yr~^, and associated stellar masses 
of the order of ~ 10^" M© per cycle. We note that the in- 
crease in mass due to the feedback-driven star formation 
is not huge, but still a non negligible fraction of the total 
stellar mass. Thus a few such episodes occurring during the 
lifetime of the galaxy could contribute to the overall stellar 
mass growth. 

In the central regions, the radiation field can be very 
strong and efficiently heat the surrounding medium via 
photoionization and Compton scattering. In particular, 
Compton heating may dominate in the innermost regions, 
where matter is highly ionized. Such a Compton-heated zone 
can extend up to a distance of ~1 kpc for a typ ical quasar 
spectral energy distribution (|Sazonov et al.ll2005l ). Thus gas 
is prevented from cooling, and consequently star formation 
is inhibited in this central region. Compton heating may 
therefore be responsible for quenching star formation in the 
inner core of the host galaxy. Indeed numerical simulations, 
including radiative feedback, show the presence of an inner 
~kpc region heated to high temperatures via Compton scat- 
tering, where the amount of cold gas is sig nificantly reduc ed 
and star formation is locally suppressed ()Kim et al.ll201ll ). 



4 SOME OBSERVATIONAL PREDICTIONS 

Here we set out a number of observational predictions of 
our model. We have seen that stars formed in the outflowing 
shell eventually follow bound orbits in the gravitational po- 
tential of the galaxy. The mean radius of the stars is mainly 
determined by the initial conditions: in general, high shell 
velocities lead to large mean radii. As a result, a significant 
increase in the size of the galaxy can be easily obtained. 
On the other hand, the associated stellar mass added in the 
process is not very large. Considering a few AGN feedback 
episodes during the lifetime of a galaxy, we would thus ex- 
pect a strong increase in size but with a modest increase 
in stellar mass. This characteristic property can be directly 
compared with the observed size and mass evolution of mas- 
sive galaxies. We should recall, however, that our model de- 
scribes active phases when the central AGN is switched on, 
while galaxies selected in observational samples are mostly 
found in the quiescent regime. 

In our picture, star formation is triggered by AGN feed- 
back and new stars are formed in the feedback-driven shell. 
We therefore expect a good correlation between star forma- 
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tion and outflow properties. We predict the star formation 
rate to scale with the outflow velocity (Eq. [3]); we also ex- 
pect more luminous, or equivalently more massive, objects 
to have higher star formation rates. The newly formed stars 
are likely to follow radial orbits, at least in the initial phases. 
Outflows are now commonly observed in galaxies, and the 
expected connection with star formation activity can be di- 
rectly searched for in the observational data. 

Star formation in the outflowing shell can be traced by 
UV and Ha emissions, characteristic of young stellar pop- 
ulations. In cases where the Ha emission can be directly 
detected, a detailed analysis of its emission line profiles and 
in particular broad wings, can give indications on the ra- 
dial motion of the stars. However, much of the intrinsic UV 
emission associated with star formation could be absorbed 
by dust grains, resulting in obscuration and reddening. We 
recall that radiation pressure on dust is the driving feed- 
back mechanism, which pushes material out to large radii, 
and thus expect significant amounts of dusty gas spread over 
the entire galaxy. We also note that the formation of new 
stars is expected to occur in dense clumps within the out- 
fiowing shell, which may travel more slowly compared to the 
mean shell speed at any radius. This suggests that dusty gas 
might travel ahead of the stars, and form a sort of screen 
shielding the intrinsic emission of young stars. We would 
therefore predict obscured star formation and dust redden- 
ing in massive growing galaxies. 

As new stars are formed at ever larger radii, a low- 
density envelope builds up in the outer regions of the host 
galaxy. The galaxy outskirts are likely to have bluer colours, 
due to the presence of young stellar populations, compared 
to the inner core, where Compton heating inhibits star for- 
mation. Negative colour gradients are then expected in the 
radial direction, with a central red core and bluer outskirts. 
But radial colour gradients might simply refiect metallicity 
gradients, and one needs to disentangle the two effects. If 
the outfiowing shell remains trapped in the galactic halo, 
gas may eventually fall back and provide fuel for further 
accretion and star formation in a subsequent episode. New 
HST observations indeed suggest that plenty of cool gas, po- 
tentially available f or re-accretion, is bound to the haloes of 
early type galaxies ijThom et al.ll20l3 'l. In a cycling scenario, 
material can be continuously recycled, leading to metal en- 
richment. Subsequent feedback cycles are fuelled by recycled 
gas which failed to undergo star formation and has fallen 
back in as well as more pristine gas from cold flows of exter- 
nal gas. This, in turn, could lead to a complex relationship 
between metallicity and feedback through the dust compo- 
nent, which is the underlying driving mechanism. 

Due to the entwinement between star formation and 
dust absorption, the detailed colour evolution of a galaxy 
is complex and difficult to predict. New stars are formed 
in the feedback-driven shell, and at the same time dusty 
gas is spread out to large radii. We thus expect intrinsically 
bluer colours in the outskirts due to the presence of young 
stellar populations, but also reddening effects due to the 
presence of dust. The dust and stellar evolutions are closely 
coupled, essentially occurring on comparable timescales, and 
the resulting colour evolution is time-dependent. 

In the above sketched picture, we have implicitly consid- 
ered the simple case of an isolated galaxy. But in a more real- 
istic situation, the galaxy is likely to reside in a group or clus- 



ter environment, where tidal interactions with neighbouring 
objects can be important. In particular, the extended low 
density envelopes of the growing galaxies can be easily af- 
fected by tidal torques, and possibly disrupted or even in- 
duced to assume some form of rotational motion by the pas- 
sage of another galaxy in the vicinity. Furthermore, if the 
outflow has a bipolar pattern, as expected if radiation is orig- 
inating from an inner accretion disc, then the growing galaxy 
might present an elongated shape, which may superflcially 
resemble a disc-like object. This may in part account for 
the disk-dominated morphologies report ed in recent obser- 
vations of high-redshift {z ~ 2) galaxies (jvan der Wei et aP 
20111), as it could be difficult to clearly distinguish the ob- 
served velocity field from rotational motion. 

Finally, since our model is based on AGN feedback, we 
expect the feedback-driven star formation to be most ef- 
fective at high redshifts, clo se to the peak epoch of AGN 
activity {z ~ 2) iMerloni fc H einz 2008). Therefore this par- 
ticular mode of star formation should mainly contribute to 
the galaxy growth in the high redshift regime (1 < z < 2). 
In this context, it is interesting to note that the peak epoch 
of both AGN and star formation activities are found to be 
roughly coincident, and it is even possible that this parallel 
evolution is triggered by AGN feedback. 



5 DISCUSSION 

In our picture, star formation is triggered by AGN feedback 
and takes place in the feedback-driven shell. We thus ex- 
pect correlations between star formation and outflow prop- 
erties. There is now ample observational e vidence of galactic- 
scale outflows detected in several sources (|Sturm et al.ll201ll : 
iTombesi et al.|[2012l . and references therein). However, con- 
trasting results concerning possible connections between 
AGN and star for mation act i vities have been reported 
in the literature JPage et al.l l2012l: ISantini et al.l I2OI2I : 
iRovilos et al] |2012| : [Harrison et al.l |2012| '1. The relation is 
usually analysed by trying to separate the dominant emis- 
sion components: X-rays are assumed to trace AGN activity, 
while far infrared/sub-millimetre emission is considered as 
a proxy for star formation. 

Based on He rschel SPIRE sub-millimetre observations, 
iPage et al.l (|2012l ) claim that star formation is suppressed 
in powerful AGNs with X-ray luminosities exceeding Lx > 
10'*"^ erg s~^ in the redshift range 1 < 2: < 3. In contrast, 
ISantini et"ai](|2012l '). using Herschel PACS observations, find 
that AGN hosts have enhanced star formation with respect 
to inactive galaxies of comparable mass in the redshift range 
0.5 < z < 2.5. Further, dividing the sample into high- and 
low- X-ray luminosities, they observe that the star forma- 
tion enhancement is more important in more luminous ob- 
jects. In a recent analysis, combining deep XMM-Newton 
and Herschel observations, a positive correlation is found 
between the specific star formation rate (sSFR) and X-ray 
luminosity for luminous objects {Lx ^ 10^"^ erg s~^) at high 
redshifts (z > 1); whereas no such correlation is detected in 
lower luminosity and lower redshift systems (|Rovilos et al.l 
2013). Another recent study, extending the analysis to a 
larger sample including the COSMOS field, show that there 
is no clear evidence of suppr ession of star formatio n in high- 
luminosity AGNs at z~l-3 ([Harrison et al.|[2012l ). Instead, 
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the mean star formation rate of luminous objects {Lx > 
lO** erg s~^) is found to be comparable to that of lower lumi- 
nosity sources, also consistent with the star formation rates 
of typical star- forming galaxies at similar redshifts ((SFR) 
-100 Mgyr"^). Although the issue is still under debate, 
there seems to be some observational evidence for a cor- 
relation between AGN and star formation activities. Such 
positive correlations are expected in the initial phases of the 
feedback-driven star formation and would be naturally ex- 
plained in our framework. 

A particular characteristic emerging from observational 
studies of the evolution of massive galaxies is the inside-out 
growth pattern. The huge masses coupled with the small 
radii observed in the high-redshift compact galaxies imply 
extremely high densities. But the central densities (mea- 
sured within a fixed physical radius of ~1 kpc) of massive 
galaxies at z ~ 2.3 are found to be quite similar to that of 
local elliptical counterparts, leading to the suggestion that 
the co mpact galaxies may fo rm the core of present-day ellip- 
ticals l|Bezanson et al.ll2009l ). Analysis of massive galaxies, 
selected at constant number density in the redshift range 
< z < 2, also show that the mass within a fixed inner 
radius of 5 kpc is almost const ant, while the main mas s 
growth takes place at larger radii (jvan Dokkum et"al]|2010l ). 
More recent results, based on Hubble Space Telescope (HST) 
imaging, confirm that most of the mass within the central 2 
kpc is already assembled at z ~ 2, and that the subsequent 
grow th occurs in the outer regions of the galaxy l|Patel et al.l 
|2012| ). All these observations seem to indicate that massive 
galaxies grow inside-out, gradually assembling their stellar 
mass around the already formed compact cores. 

Another important feature is that the evolution in size 
and mass are not parallel: the increase in size is much 
stronger than the increase in mass. The relation between the 
effective radius and the stellar mass can be parametrized 
as Re oc M" with a > 2 (jPa tel et al. 2 012,). In addition, 
the b ulk of the mass and size evolution seem to occur at 
z>l l|Perez-Gonzalez et al.ir2008l : lFan et al.ll2010l'l. with the 
grow th rate being faster at higher redshifts ( Newman et al.l 
|2012| ). 

As mentioned in the Introduction, different physical 
mechanisms have been proposed in order to explain the ob- 
served size and mass evolution of massive galaxies. In the 
puffing-up scenario, quasar feedback removes large amounts 
of gas from the central regions, l eading to an expansion of 
the system at fixed stellar mass (|Fan et al.ll2008 . 2010). As 
there is no increase in mass, while the size increases signif- 
icantly, the velocity dispersion of the galaxy is expected to 
decrease considerably. In contrast, observational measure- 
ment s indicate only a mild ev olution in the velocity disper- 
sion (|Cenarro fc Truiilld 120091 . The puffing-up scenario is 
therefore not supported by observations. 

In major mergers of two equal mass galaxies, the con- 
servation of kinetic energy implies Re (x M, and the velocity 
dispersion remains constant. This indicates a linear increase 
in both size and mass, e.g. a doubling in size requires a dou- 
bling in mass. If the observed size growth is attributed to 
major mergers, the corresponding mass growth would yield 
too many massive galaxies (~ IQi^^Mq), in contradiction 
with the number density observed in the prese nt day uni- 
verse (IBezanson et al. I l2009l : iHopkins et al]|2010l ) . Moreover, 
direct estimates of the pair fraction of massive galaxies sug- 



gest that the merger fraction is not sufficient to account for 
the o bserved size growth from z ~ 2 to z ~ (|Man et al.l 
I2012D . Therefore major merging can be ruled out as the main 
growth mechanism for massive galaxies. 

In the case of minor mergers, the primary galaxy ac- 
cretes material by merging with lower mass companions 
(typically with ratio 1:10). Conservation of energy implies 
that the effective radius is proportional to the square of 
the mass. Re oc . Thus for a given increase in mass, 
the corresponding increase in size is much more efficient 
compared to the case of major mergers. Numerical simu- 
lations of minor mergers also show that the accreted low- 
density material tends to settle in the outer parts of the 
primary galaxy, implying a gr owth pattern consistent with 
the in side-out growth scenario l|Naab et al.ll2009l :[o ser et al.l 
l2012h . Therefore the minor merger model is currently con- 
sidered as the favoured interpretation in accounti ng for 
the observed evolution of massive galaxies (Bczans on et al.l 
I2OO9I : iHopkins et alll2010l ) . However, recent results from the 
CANDELS survey suggest that minor mergers may only ex- 
plain the size growth at z < 1; whereas the faster growth 
rate observed at higher redshifts cannot be a ccounted for, 
and requires additional physical processes ((Newman et al.l 
l2012h . 

In our model, new stars are formed in the AGN 
feedback- driven shell at increasingly larger radii and pop- 
ulate the outer regions of the host galaxy. It is also possible 
that some of the stars will be pushed beyond the gravita- 
tional potential of the host and be ejected from the galaxy. 
Such young stars distributed in the halo may form part of 
the intergalactic stellar population and contribute to the en- 
richment of the intracluster medium. We have seen that the 
increase in galaxy size is accompanied by a modest increase 
in stellar mass, in agreement with the observational trend. 
The succession of several such feedback-driven star forma- 
tion episodes may thus account for the evolution of massive 
galaxies. Most of the galaxy growth is observed to occur at 
large radii, while the c e ntral region remains roughly constant 
jBezansonet al.ll2009l : Ivan Dokkum et akllioiol : iPatel et al] 
l2012ll . The absence of significant star formation in the cen- 
tral core, due to Compton heating, may explain the similar- 
ity of the core component observed in massive galaxies. 

Massive, compact galaxies at z ~ 2 are often observed 
in a quiescent phase. Deep spectroscopic observations of a 
typical compact quiescent galaxy at z ~ 2.2 suggest that 
quenchin g of star formation can be very effective at these 
redshifts (|Kriek et all 120091 '). As a result, a red-sequence is 
likely to be already in place at z ~ 2, with galaxies mainly 
in a post-starburst phase. However, the colour evolution of 
the red-sequence to the present cannot be explained by just 
assuming passive evolution: simple aging of the stellar pop- 
ulations would produce galajcies at z ~ that are redder 
than observed (|Kriek et al.ll200^ ). New starburst events oc- 
curring in red-sequence galaxies at later times might help in 
reducing the strong colour evolution. In our picture, several 
episodes of star formation associated with the AGN feedback 
activity are envisaged over the lifetime of a galaxy. As young 
stars are formed at increasingly larger radii in the outflow- 
ing shell, we also expect a radial colour gradient to develop 
within the galaxy, leading to bluer outskirts. Indeed negative 
colour gradients are observed in massive early type galaxies 
in the redshift range 1.3 < z < 2.5, with colours becoming 



© 2012 RAS, MNRAS 000, [THi] 



6 W. Ishibashi, A. C. Fabian and R. E. A. Canning 



bluer with increasing radius (|Guo et al.ll2oTll ). Similar neg- 
ative colour gradients in the radial direction are detected in 
a sample of galaxies a t z ~ 2 with a variety of morphologies 
l|Szomoru et al.ll201ll ). 

Compact galaxies at high redshifts are observed to be 
'red and dead', with little gas and no ongoing star forma- 
tion. It is plausible that an earlier episode of AGN feedback 
has removed cold gas from the core of the galaxy, but gas 
should still be present in the outskirts of the galaxy and in 
the surrounding circumgalactic medium, trapped in the dark 
halo. Indeed large reservoirs of cool gas are observed in th e 
haloes of early-type galaxies at z < 1 (|Thom et al.ll2012l ). 
and recently a significant amount of cool gas extending to 
hundreds of kilopars ecs has also been dete cted in quasar 
host haloes at 2 ~ 2 (|Prochaska et al]|2013l ). Over the cos- 
mic time, from z ^ 2 to the present, we expect galaxies to 
undergo several gas re-accretion phases, which should power 
AGN activity and associated feedback-driven star formation 
cycles. The observed reservoirs of gas surrounding the pro- 
genitors of today's ellipticals may provide the required gas 
fuel. 

A significant amount of dust [Ay > 1) is also required 
in order to efficiently drive AGN feedback. Recent results 
based on Herschel observations of the Crab Nebula suggest 
that the producti on of dust in core-c ollapse supernovae can 
be quite efficient (| Gomez et al.|[20l3 ). This implies that new 
dust can be produced by supernova explosions and spread 
into the surrounding star-forming medium, enhancing the 
overall feedback process. At the same time, dust may also 
be responsible for extinction and reddening. We thus expect 
dust-obscured star formation in the host galaxy. The ob- 
served red colours of massive galaxies at high redshifts are 
usually interpreted as due to old stellar populations, with 
negligible star formation. But dust reddening may also play 
a non negligible role, and the observed red colours might 
be partially attributed to dust effects. In fact, samples of 
galaxies at 1.5 < 2 < 3 selected from the GOODS-NICMOS 
survey show that massive galaxies have high star formation 
rates, of the order of ~100 M(7,yr~ ^ , but are still dom inated 
by red rest-frame (U-B) colours l|Bauer et al.l I2OI1I '). This 
suggests that star formation in most of these systems is 
obscured by dust. A number of heavily reddened quasars, 
requiring large dust extinc tion (Ay ^ 3) have also been re- 
cently observed at 2 ~ 2 (|Banerii et al.|[201^ '). Thus there 
seems to be some observational evidence for dusty star for- 
mation and dust reddening occurring in massive galaxies at 
high redshifts. 

Different forms of star formation triggering have been 
previously discussed in other contexts. Early studies sug- 
gested triggered star formation, associated with the propa- 
gation of relativistic jets, as a possible explanation for the 
observed alignment of optical and radio emissions in radio 
galaxies (Rocs 198^). A jet-triggered positive feedback can 
also enhance the star forma tion in luminous starbursts (jSilkl 
120051 ). ISilk fc Nusseij l|2010l ') invoke star formation triggered 
by AGN winds to boost feedback, from the combination of 
both AGN and supernovae, in order to drive gas out of the 
host galaxy. More recently, detailed numerical simulations 
of radio mode feedback, with jet-induced sta r formation, 
have also been performed l|Gaibler et al.l [20121 ). While jet- 
triggered star formation is mainly associated with radio-loud 



objects, the radiation pressure- driven feedback considered 
here may be applied to a wider range of sources. 

The particular galaxy growth mechanism discussed here 
mainly operates at high redshifts, where AGN activity is 
more pronounced, and where the rate of minor mergers 
seem to be insufficient. Therefore AGN feedback-driven star 
formation may explain the size growth at high redshifts 
(1 J; ^ 2), while minor merging may play a major role 
at lower redshifts (0 < 2 < 1), with a combination of differ- 
ent physical mechanisms dominating at different epochs. 
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